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Carl J. Daniele and Fred Teren 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 


Abstract 

A simulation technique is presented for the 
prediction of compressor stall for axial-flow com- 
pressors for clean and distorted inlet flow. The 
simulation is implemented on the digital computer 
and uses stage stacking and lumped- volume gas dy- 
namics. The resulting nonlinear differential 
equations are linearized about a steady-state op- 
erating point, and a Routh-Hurwitz stability test 
is performed on the linear system matrix. Par- 
allel compressor theory is utilized to extend the 
technique to the distorted inlet flow problem. 

The method is applied to the eight-stage J85-13 
compressor. Analytical stall prediction for the 
undistorted stall line shows good agreement with 
experimental results. The predicted stall line 
for distorted inlet flow is in agreement with ex- 
perimental results only for large distortion ex- 
tents and/or low distortion levels. Results for 
low distortion extents and high distortion level 
do not agree with experimental results. 

Nomenclature 

S& area, m 2 ; ft 2 

C coefficient, rpnT/R 2 ; rp m 2 /°R 2 

CONST constant for weight flow derivatives, 

N/m 2 ; lbf/ft 2 

g gravitational constant, 1. O(kg-rn)/ (N-sec 2 ); 

32 . 17(lbm-f t)/ (lbf-sec 2 ) 

K coefficient, rpm-sec/m; rpm-sec/ft 

KCHOK choked flow coefficient, (kg-/K)/ (N-sec) ; 
(lbm-/°R) / (lbf -sec) 

L length, m; ft 

N rotational speed, rpm 

P pressure, N/m 2 ; lbf/ft 2 

P2 input pressure, N/m 2 ; lbf/ft 2 

R universal gas constant, 287 (N-m) / (kg-K ) \ 

53 . 3(lb£-f t) / (lbm-°R) 

T temperature, K; °R 

T2 input temperature, K; °R N 

t time , sec 

V volume, m^; 

v velocity, m/sec; ft /sec 

W weight flow, kg /sec; lbm/sec 

v ratio of specific heats 


6 ratio of total pressure to sea level 

pressure 

0 ratio of total temperature to sea level 

temperature 

p weight density, kg/m^; lbm/ft^ 

$ flow coefficient 

pressure coefficient 
temperature coefficient 

Subscripts : 

b variable associated with stage bleed 

c variable associated with stage charac- 

teristic 

h variable associated with undistorted 

side of compressor 

A variable associated with distorted side 

of compressor 

n stage number designation 

OGV outlet guide vanes 

s static condition 

t total condition 

v variable associated with stage volume 

z axial direction 

1,2,..., 8 stage numbers 

9 burner volume 

10 nozzle exit 
Superscripts : 

derivative with respect to time 

maximum value of quantity used for 
normalizing 

I. Introduction 

The transient performance of airplane propul- 
sion systems using axial-flow compressors is 
greatly affected by the compressor stall charac- 
teristics. Compressor stall adversely affects 
engine performance by reducing thrust and engine 
airflow, and can cause inlet unstarts in supersonic 
aircraft. The stall characteristics of axial-flow 
compressors are directly affected by inlet distor- 
tion. Thus, the need for determining the stall 
line for these compressors under both distorted and 
undistorted inlet flow conditions is evident. 



Experimental methods are generally used to 
determine stall lines for axial-flow compressors. 
However, engine experiments to define stall lines 
are costly and time consuming; thus, mathematical 
techniques to predict compressor stall conditions 
are being sought and developed. The procedures in 
use at present are one -dimensional and require well 
defined stage pressure and temperature character- 
istics. These characteristics must be determined 
experimentally. One method used to predict com- 
pressor stall is presented in reference 1, in which, 
the undistorted stall line for the J85-13 axial- 
flow compressor was determined from an analog simu- 
lation and compared to experimental results. Re- 
cently , this work has been extended to digital com- 
puters. Both simulations employ compressor stage 
stacking and lumped volume gas dynamics. On the 
analog computer, a limit cycle or instability of 
the time varying solution is correlated with ex- 
perimentally observed compressor stall. The digi- 
tal method involves linearization of the nonlinear 
system of equations (state equations) which defines 
the compressor simulation about a steady-state op- 
erating point, thereby obtaining a linear system 
matrix. System stability at this compressor oper- 
ating point is then determined by a Routh-Hurwitz 
test of the system matrix, using digital routines 
found in reference 3. In reference 2, the stall 
line for the J85-13 was determined by the digital 
stability analysis. Comparison with the stall line 
generated by the analog simulation from reference 1 
showed good agreement. 


as in reference 2. With this model, the predicted 
compressor stall point depends on the operating 
point of the undistorted side as well as the dis- 
torted side. 

Results are presented first which compare the 
analytically predicted stall line with experimental 
data for undistorted flow. Then using a parallel 
compressor model and the digital stability analysis, 
stall points are predicted for varying distortion 
extents and intensities and are compared with ex- 
perimental results. 

II. Engine Description 

The analysis in this report is applied to the 
J85-13 afterburning turbojet engine. The engine 
consists of an eight-stage, axial-flow compressor, 
a two-stage turbine, variable inlet and outlet guide 
vanes, controlled interstage bleeds, and variable 
exhaust nozzle. The inlet guide vanes are scheduled 
as a function of corrected speed; the outlet guide 
vanes are fixed, but the pressure drop across them 
varies with corrected weight flow. The Interstage 
bleeds are on the third, fourth, and fifth com- 
pressor stages and are also scheduled as a function 
of corrected speed. (The bleeds are closed above 
94.0 percent corrected speed, are full open at 
80 percent speed and below, and vary with speed in 
between, ) 

III. Simulation Model: Undistorted Flow 


A technique which is sometimes useful in pre- 
dicting stall for distorted inlet flow is parallel 
compressor theory (first postulated in ref. 4). 

The theory assumes that a compressor with distorted 
inlet flow can be represented by two fictional 
compressors, one having an inlet pressure equal to 
the distorted or low total pressure, and the other 
having an inlet total pressure equal to the undis- 
torted or high total pressure. The distortion at 
the inlet is assumed to have a square wave form. 
Both compressors are assumed to discharge to a 
common static pressure and operate on the previ- 
ously determined clean inlet speed line. No cir- 
cumferential flow is exchanged between the two 
compressors. Stall is assumed to occur when the 
pressure ratio of the "low compressor 11 is equal to 
the clean inlet stall pressure ratio. Prior knowl- 
edge of an experimentally determined clean inlet 
stall pressure ratio is therefore necessary. 

This theory has been investigated in refer- 
ences 5 through 8. Xn these references it is ob- 
served that parallel compressor theory gives accu- 
rate stall prediction for large distortion extent, 
but becomes progressively more conservative as ex- 
tent decreases. 


The single compressor with undistorted inlet 
flow is modeled as follows. The multistage com- 
pressor portion of the simulation uses steady-state 
stage characteristics and an associated stage vol- 
ume for gas dynamics to represent individual stages 
Each compressor stage is modeled with pressure and 
temperature maps; each volume, with lumped conti- 
nuity, energy, and momentum equations. The inlet 
and outlet guide vanes are modeled as pressure 
drops into and out of the first and last stages , 
respectively, the amount of pressure drop being a 
function of corrected speed. Bleed flows are in- 
cluded for the third, fourth, and fifth stages; the 
amount of bleed flow is also a function of cor- 
rected speed. The bleeds are assumed choked. Be- 
tween the last compressor stage and turbine, a 
burner volume is modeled. No heat addition is in- 
cluded, however, because only the compressor, 
rather than the whole engine, is of interest for 
this report. Behind the burner a choked nozzle is 
used to terminate the simulation. The choked noz- 
zle simulates a choked turbine in the actual engine 
A schematic of an ideal compressor stage is shown 
in figure 1. 

Steady-State Equations 


In this paper a method is described for pre- 
dicting stall line location for circumferentially 
distorted inlet flow conditions which combines the 
digital stability analysis of reference 2 with a 
parallel compressor model. The objectives of this 
analysis are to eliminate the need for a knowledge 
of a reference clean inlet stall pressure ratio 
and to improve the accuracy of stall line predic- 
tion. The digital simulation as in reference 2, 
used stage stacking and lumped volume gas dynamics. 
Two separate compressors were simulated with a 
static pressure balance in the burner volume at 
steady- state. The combined nonlinear system of 
equations is linearized and checked for stability 


The steady-state performance of the axial-flow 
compressor is calculated by using the stage stack- 
ing technique described in reference 1. In this 
technique the discharge conditions of one stage are 
used as the inlet conditions to the next stage. 

The steady-state performance for each stage is 
given by pressure and temperature rise coefficients 
and These coefficients are defined as 
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with the temperature in the burner volume assumed 
(2) equal to the temperature In the outlet guide vanes: 


Represent at ive maps of ^ and ^ as a function 
of flow coefficient $ are shown in figure 2, The 
flow coefficient is defined as 
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Thus, the undistorted compressor simulation con- 
sists of 26 state equations - three for each of the 
eight compressor stages and two for the burner 
volume . 


The steady -st ate performance of the burner volume 
is modeled by assuming that the pressure, tempera- 
ture , and weight flow into and out of the volume 
are the same. The nozzle and bleed ports are as- 
sumed choked, with the amount of bleed flow area a 
function of corrected speed. The choked flow equa- 
tion used is 


W = KCHOKj-r — - (4) 


IV. Simulation Model: Distorted Flow 

Parallel compressor theory is used as a basis 
for extending the undistorted inlet flow model to 
distorted inlet flow conditions. A schematic dia- 
gram of the parallel compressor model is shown in 
figure 3. The assumptions made are that two sepa- 
rate compressors are modeled with different total 
pressures at their inlet and that no crossflow 
occurs between the compressor segments. At steady 
state, static pressure is balanced at the discharge 
of the compressor segments. 


The steady-state equations are used to determine an 
operating point on a speed line. 

Dynamic Equations 


Compressor dynamics are simulated by using a 
lumped volume for each stage. Each lumped volume 
is located just downstream of the stage’s active 
element. The volume dynamics are modeled through 
the application of continuity, energy, and momen- 
tum equations. The equations are 


_d_ 

dt 


sv.n' 


- W 

c^n+l' 


( 5 ) 


Steady-State Equations 

Two completely separate compressors are 
modeled with different total pressures at their in- 
lets. The steady-state performance of each com- 
pressor is calculated using equations (1) to (3) 
for the eight compressor stages and the stage maps 
determined from clean inlet flow experimental tests 

However, additional equations are added in the 
burner volume to account for the different exit 
conditions for each compressor segment. In the 
burner volume, the weight flows from the two com- 
pressors are added and the exit temperatures are 
mass averaged, thus 


dt ^ P sv,n^tv,n^ 


( 12 ) 


F [T tc,n W c,n T tv,n (W c,n+l + “b.n^ 


( 6 ) 


and 


T _ (M c,<H T tv,9* + W e,9h T tv,9h^ fl _. 

*v,9 * 7* 

c,9H + c,9h 




dt " — (P tc,a ' P tv,«> 


(7) 


While the equation of state for the gas in the 
volume is 


The total pressure in the burner volume is calcu- 
lated by using the static pressure and the compres- 
sible flow equations. Once the pressure, tempera- 
ture, and weight flow are calculated, a choked noz- 
zle equation (4) is used to terminate the model. 

Dynamic Equations 


P = Rc T 
tv,n ^sv,n tv,n 


( 8 ) 


No dynamic effects are modeled for the inlet and 
outlet guide vanes, interstage bleeds, or choked 
nozzle . A lumped volume is included In the burner. 

The dynamic equations used to model the burner 
volume are 


For both compressors simulated, the dynamic 
equations for the single compressor (eqs. (5) 
to (8)) apply for the eight compressor stages. In 
these equations, however, the volumes and areas are 
sized to account for the distortion extent. Equa- 
tions (9) and (10) used to simulate dynamics in the 
burner volume are modified to account for the par- 
allel compressors by 


dt (p sv,9 T tv,9 ) “ V,/ l ”c,9 ”c,10 ] 


where Is the choked nozzle flow, and 


a . 8 

dt (W c,9> =L7 (P tc,9 - P tv,9> 


( 10 ) 


dt (p Bv,9 T tv,9 ) = T ty,9 ( ”c,9«. + ”c,9h ~ 


(14) 


where Win is the choked nozzle flow. Also 


dt (W c,9h> 
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and 


£ » tc ,9l - P tv,9 - °»»V < 16 > 

where CONST is a steady-state pressure drop and 
is set such that d/dt is equal to zero in 

steady state. The temperature in the burner volume 
is calculated from equation (13) . Thus , the par- 
allel compressor model consists of 51 state equa- 
tions - twenty four for each compressor segment and 
three for the burner volume. 

V, Stability Analysis 

The most direct way to determine the stability 
of a system of nonlinear differential equations is 
to obtain the transient response to a small dis- 
turbance, starting at a steady-state operating 
point. A stable system will return to the steady- 
state condition once the disturbance has subsided. 
An unstable system, however, will either diverge 
or undergo limit cycle oscillations. This method 
of stall prediction was used in reference 1 wherein 
an analog computer simulation was used to obtain 
the transient solution and predict compressor stall 
for clean inlet flow for the J85 engine. 

It was desired to extend the stability anal- 
ysis technique of reference 1 to the prediction of 
compressor stall for distorted inlet flow. This 
requires the simulation of two separate compres- 
sors, as discussed earlier. It was not possible to 
construct such a simulation on the analog computers 
at the Lewis Research Center, because the large 
amount of equipment required exceeded the amount 
of available analog equipment. Therefore, it was 
decided to construct an equivalent simulation on 
the digital computer. Results from this simula- 
tion (unpublished) showed that a time step of 
50 microseconds or less was required to avoid 
numerically induced instability. As a result, 
about one hour of computer time was required for 
1 second of real-time data. Because of this, the 
computer time required to conduct substantial sta- 
bility analysis on the digital computer was pro- 
hibitive. 

In order to make use of the digital computer 
for stability prediction and avoid large computer 
times, a new stability analysis technique was de- 
veloped and is the subject of reference 2. This 
technique does not require numerical integration of 
differential equations. In addition to stability 
prediction, the method is very useful and efficient 
for the calculation of frequency responses, eigen- 
values, and eigenvectors, as shown in reference 2. 

The stability analysis in reference 2 is ap- 
plied to the model by first determining a steady- 
state point on a speed line. The nonlinear system 
of equations (state equations) which defines the 
compressor simulation about the operating point is 
linearized and a linear system matrix is obtained. 
The elements of the system matrix are the partial 
derivatives of the state equations with respect to 
the state variables. The linearization of the 
state equations and calculation of the matrix ele- 
ments are accomplished by finite differences on the 
digital computer. The stability of the system 
matrix is checked by a Routh-Hurwi t z test. When 
this test indicates Instability, the linear system 
is unstable and the nonlinear system is also un- 


stable through application of the first method of 
Lyapunov.^ 


The stability analysis is first applied to the 
undistorted single compressor model. A point is 
chosen on a speed line and its stability determined 
by the Routh-Hurwi tz test. If it is stable, points 
are chosen at successively lower flow rates on the 
speed line until instability is indicated. The sta 
bility analysis is also applied to the prediction 
of stall for distorted inlet flow by using the 
51-variable dynamic model described earlier. This 
method should improve on classical parallel com- 
pressor theory since it does not require an experi- 
mentally determined clean inlet stall point as ref- 
erence. Also there should be an improvement in 
overall stall prediction since the whole compressor 
model (both imdistorted and distorted segments) is 
considered in determining compressor stability. 


V. Results 


Results are presented first which compare the 
analytically predicted stall line with experimental 
data for undistorted flow. Then using a parallel 
compressor model and the digital stability analysis 
stall points are predicted for varying distortion 
extents and intensities and compared with experi- 
mental res \ilts. 

Data Base 


The experimental data used is from a recent 
test of a J85 engine in a Lewis altitude facility 
where stall line definition for both undistorted 
and distorted inlet flow was accomplished. Stage 
temperature and pressure maps were determined and 
used in the simulation. These stage maps give a 
mean line representation of the stage 
characteristics - data is taken at the stator exit 
at three radial positions and those results are 
area weighted radially to get mean line values. 

From these values one set of pressure and tempera- 
ture stage maps is determined. This set of maps is 
used for all compressor speeds. 

Undistorted Stall Line 


The overall compressor map as determined by 
the simulation Is compared to experimental results 
in figure 4. This figure shows overall compressor 
ratio as a function of weight flow corrected to the 
compressor face. Speeds of 80, 87, 94, and 
100 percent are shown. The circles represent ex- 
perimental data points and the solid lines are the 
analytical speed lines. Darkened triangles and 
circles represent the analytical and experimental 
stall points, respectively. The dashed line is the 
experimental stall line. Stall points predicted by 
the stability analysis agree well with experimental 
results. The analytical speed lines agree well 
with the experimental speed lines except at 94 per- 
cent. However, this speed line differs from pre- 
vious J85 experimental data and is suspect. The 
fact that the interstage bleeds begin to open at 
94 percent speed may have some effect on the shape 
of the speed line and warrants further investiga- 
tion. Therefore, in this report, further results 
for the 94 percent speed line will be included but 
not discussed. 

Distortion Results 


The ejcperimental data base for the analysis 
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consisted of four different density screens and 
four different distortion extents. Table 1 lists 
the different combinations of screens and extents 
used. Figure 5 shows the ability of the analysis 
to predict performance and stall limits. Fig- 
ure 5(a) is for the ISO degree extent and 5(b) for 
the 30 degree extent. In both cases the results 
shown are for the medium density screen. Good per- 
formance and stall prediction is shown for the 80 
and 100 percent speed lines in figure 5(a), Stall 
prediction for the 87 percent speed line is fair 
but there is good agreement with the experimental 
speed line . 

In figure 5(b) there is good agreement between 
analytical and experimental speed lines. However, 
stall prediction is conservative in all cases with 
the 100 percent speed line being more acceptable. 
The analytical speed lines have also been calcu- 
lated beyond the stall points, but the results are 
not shown here. The analytical speed lines diverge 
from the experimental results for weight flows be- 
low the predicted stall point. 

Figure 6 shows the effect of varying distor- 
tion extent with constant screen density. Only 
experimental and analytical stall points, and not 
speed lines, are shown. Open symbols are analyti- 
cal results and darkened symbols are experimental 
results. Figures 6(a), (b), and (c) show results 
for the high, medium, and low density screens, re- 
spectively. Figure 6(a) shows good agreement for 
the 180 degree extent but the agreement deterio- 
rates at 90 and 60 degrees. As the screen density 
decreases (figs. 6(b) and (c)) the agreement for 
the small extents becomes better. Also, from fig- 
ure 6(a), experimental results indicate (neglecting 
scatter) that the stall limit increases as extent 
decreases fox fixed screen density as expected. 
However, the analytical results indicate the oppo- 
site effect. 

Figure 7 shows the same data as figure 6 but 
is plotted to show the effect of screen density as 
distortion extent is held constant. In this fig- 
ure also, only the experimental and analytical 
stall points are shown. Figure 7(a) shows that for 
the 180 degree extent , there is good agreement be- 
tween the analytical and experimental results. 
Figures 7(b) and (c) with 90 and 60 degree extents, 
respectively, show increasing degradations in stall 
prediction. In all cases, however, the expected 
trend of decreasing compressor performance with 
increasing distortion level is present in the ana- 
lytical results and also in the experimental re- 
sults. 

VII. Concluding Remarks 

Presented in this report is the application of 
the compressor stall prediction method from ref- 
erence 2 to an axial-flow compressor recently 
tested in a Lewis altitude facility. The test in- 
volved stall line definition for the compressor 
with both undistorted and distorted inlet flows. 

The stall prediction correlates well with experi- 
mental results for clean inlet flow. The inclu- 
sion of the stability analysis with the parallel 
compressor model used to simulate distorted inlet 
flow, however, produces mixed results. The corre- 
lation of instability with experimental stall is 
acceptable at 180 degree extent with all distor- 
tion levels and at smaller extents for small dis- 
tortion levels. The correlation is not good for 


extents less than 180 degrees for large distortion 
levels. This is the same disadvantage as noted in 
references 5,6, and 7 for stall determination 
using pure parallel compressor theory. Thus, no 
Improvement has been made in compressor stall pre- 
diction with the addition of the stability analysis 
to the parallel compressor model. 

The failure of the method to accurately pre- 
dict stall for distortion extents less than 180 de- 
grees may be related to the lack of coupling of the 
distorted and undistorted sectors in the compressor 
section itself. That is, as the parallel compres- 
sor model is presently incorporated in the analysis, 
the undistorted and distorted sectors both contrib- 
ute to the stall prediction method, but only through 
the burner volume. Since the burner volume is much 
larger than the stage volumes, its dynamic contri- 
bution relative to the whole system is small. 
Therefore, one obvious modification to the analysis 
would be to permit interaction between undistorted 
and distorted sectors at each stage through circum- 
ferential crossflows in the stage volumes. Thus, 
the ability of the undistorted sector to provide a 
stabilizing influence on the distorted sector would 
be enhanced. The authors are currently modifying 
the simulation described in this report to allow 
for circumferential crossflows. 
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Table I Combinations of screen densities 


and distortion extents 


■ 

Screen 

density 

. (p2 h - P V 

* 

80% speed 




% P 

average 

180 


D 

1 

Low 

2.083 


X 

X 



Medium 

4.234 


X 

X 

X 

X 

High 

4.741 


X 

X 

X 


Highest 

7.163 




X 



X denotes experimental data taken. 


The distortion levels given are indicative of the distortion 
levels for the screens. Slight differences in these numbers 
result at different speeds and along each speed line. 
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Figure 2. - Representative 
compressor stage pres- 
sure and temperature 
maps. 



Figure 3. - Parallel compressor model. 
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Figure 4. - Comparison of analytical and experi- 
mental stall lines. Undistorted flow. 
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Figure 5. - Comparison of experimental and analytical 
results. Distorted inlet flow. 
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(C) LOW DENSITY SCREEN. 


Figure 6. - Comparison of analytical and experimental stall 
points. Screen porosity constant. 
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Figure 7. - Comparison of analytical and experimental stall 
points. Distortion extent constant. 
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